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We study the step meandering instabihty on a surface characterized by the ahernation of terraces 
with different properties, as in the case of Si(OOl). The interplay between diffusion anisotropy and 
step stiffness induces a finite wavelength instability corresponding to a meandering mode. The insta- 
bility sets in beyond a threshold value which depends on the relative magnitudes of the destabilizing 
flux and the stabilizing stiffness difference. The meander dynamics is governed by the conserved 
Kuramoto-Sivashinsky equation, which display spatiotemporal coarsening. 
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Molecular beam epitaxy (MBE) is often used to grow 
nanostructures on vicinal surfaces of semiconductor and 
metallic crystals Q, |^ |^ ^ |^ . Under nonequilibrium 
growth a very rich variety of crystal surface morphologies 
are experimentally observed resulting from the nonlinear 
evolution of step bunching and meandering instabilities 

0, • Self-organized patterns arising from these insta- 
bilities may be used for the development of technological 
ap plic ations such as quantum dots and quantum wells 

The step meandering instability was originally pre- 
dicted theoretically by Bales and Zangwill jllj for a vic- 
inal surface under growth. Its origin comes from the 
asymmetry between the descending and ascending cur- 
rents of adatoms. Nonlinear extensions of this work 
have shown that the meander evolution can be de- 
scribed by amplitude equations showing diverse behav- 
iors: spatiotemporal chaos governed by the Kuramoto- 
Sivashinsky equation in the case of the Erlich-Schwoebel 
effect with desorption [l^ : nonlinear coarsening in the 
case of negligible desorption 0, ; and interrupted 
coarsening when two-dimensional anisotropy is included 
Step meandering on a Si(OOl) vicinal surface 
can also be driven by a drift electromigration current in 
the presence of alternating diffusion coefficients, even for 
symmetric adatom attachment to the steps 0, 0| . 

Recent experiments [l^, 0, on the growth 

of Si(OOl) have revealed the existence of a step bunch- 
ing instability and the development of a transverse two- 
dimensional complex structure (ripples), possibly rem- 
iniscent of a meandering instability. We have recently 
shown that the observed step bunching instability is due 
to the interplay between the elastic interactions and the 
alternation of the step parameters [j^ . In order to un- 
derstand the roughening of the Si(OOl) surface during 
growth it would be important to know if step bunch- 
ing and meandering can arise simultaneously. So far, no 
conclusive experimental evidence for the presence of the 
Schwoebel effect leading to meandering instability has 
been established. 
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FIG. 1: Sketch of the Si(OOl) vicinal surface showing the al- 
ternation of terraces and steps Sa and S^- Lines on terraces 
indicate the privileged diffusion directions. Da, Db and 7a, 
7b are the surface diffusion and step stiffness coefBcients, re- 
spectively; x", Xb and ^a{y,t), ^biVji) are step positions and 
the corresponding perturbations. 



In this Letter, we show that the difference in step stiff- 
ness and diffusion anisotropy induce a meandering insta- 
bility which can account for the complex step structures 
observed in Si(OOl) epitaxial growth. We first present 
the linear two-dimensional stability analysis of a train of 
steps using a simplified extension of the model studied 
in reference We show that the interplay between 

diffusion anisotropy and the step stiffness effect under 
growth conditions induces a finite wavelength instability 
which is maximum for the in-phase modes. We present 
the stability phase diagram in the parameters space. Our 
results are complemented by a weak nonlinear analysis of 
the step meander which reveals a coarsening dynamics. 
Using a simple similarity argument, we show that the 
characteristic coarsening exponent is 1/2 and that the 
general solution of the CKS equations can be thought as 
a superposition of parabolas. Finally we conclude this 
Letter by discussing the possible relevance of our model 
to the experiment on Si growth and the consequences of 
the simultaneous existence of bunching and meandering 
instabilities. 
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The Si(OOl) vicinal surface consists of a periodic se- 
quence of terraces where rows of 2 x 1 dimerized adatoms 
(terrace of type a) ahernate with 1x2 dimerized adatoms 
(terrace of type 6), as shown in Fig. ^ where we also in- 
troduce several notations. On the reconstructed surface 
adatoms diffuse preferentially along dimer rows, giving 
rise to an anisotropic diffusion. Therefore, the steps sep- 
arating the terraces are of two kinds. The Sa steps are 
rather straight while the Sb ones are very corrugated [i^ . 
We shall allow each step to have a different step stiffness 
coefficient 7a and 7b |26l |. For simplicity we neglect elas- 
tic interactions between steps and assume that the des- 
orption of adatoms is negligible; we also neglect Erlich- 
Schwoebel effects. Let us denote by a;"(y,i) and x]!;{y,t) 
the positions at time t of steps Sa and St respectively 
(cf. Fig. During growth, the adatom concentrations 
on each terrace C^{x,y,t) and CJ^{x,y,t), obey the fol- 
lowing diffusion equations |2^ : 



D.df.C:! = -F, 



Dad'yCl' = -F, 



(1) 
(2) 



where Da and Di, are the diffusion coefficients, and F the 
deposition flux. C" and are the difference of adatoms 
concentrations with respect to the uniform equilibrium 
concentration Co (taken to be the same on both terrace 
types). We assume for and the following bound- 
ary conditions: 
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(3) 
(4) 



which correspond to instantaneous attachment kinetics, 
and can be considered as the simplest ones capturing the 
main physical effects (no diffusion along steps and negli- 
gible transparency). The adatom equilibrium concentra- 
tions Cr^ „ and f, depend on the step curvatures 



and k]^ 11 



(5) 



with Ffc = fl^b/kBT and F^ = fl'ja/kBT (Q is the unit 
atomic surface, T the temperature, and ks the Boltz- 
mann constant). The normal velocities of each step are 
< = ia/(l + idyx^af)^/^ and < = i^/(l + {dyX^)Y^\ 
where 

x: [{Dbd^c^^ - Da{dyxi)dyC^)- 

{Dad^C: - Db{dyXl)dyC:%^^^^ , (6) 

±1 =n [{Dad,c:+' ~ Db{dyxi)dyc:+')- 

(Db - Daidyxl)dyC^)l^^,. . (7) 

In order to get a nondimcnsional version of equations 
CQ)-©, we set the unit of length to be the initial size 
of the terrace Iq and the unit of time l^/ (CoTa^Da)- 
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FIG. 2: Stability diagram in the plane (/o,5o), for cto > 1- 
The thick gray line fo — foc given by Eq. 0, separates the 
unstable region (right side), from the stable one (left side). 
The dispersion relation a{q) with its two branches is shown 
in each region. 



The system is controlled by three independent positive 
nondimcnsional parameters: 



/o = 



CoTaDa 



F>b c - Ff, 

"0 = 75- , Oo = — , (8) 



which respectively relate to the flux, diffusion anisotropy 
and step stiffness difference. 

We investigate now the linear stability of a train of 
equidistant steps traveling at a constant velocity /o when 
perturbed transversally. The shape of the steps can be 
decomposed in Fourier modes of the form x"{y,t) = 
fot + 2n + e (y , t) and x^{y,t) ^ fot + 2n + 1 + 4" (y , t) 
with and the perturbation amplitudes varying as 
exp(cr(g, (/))t -|- iqy + in(f>), where q is the wavenumbcr 
and (p the phase (see Fig. Inserting these expres- 
sions into (E3)-(I7|), we obtain the general dispersion re- 
lation a — (T(q, 0). The dispersion relation possesses, 
for each (f), two branches corresponding to a stable as 
and an unstable mode. The maximum growth rate is 
reached for the in-phase perturbation (j> = (Fig. [21. In 
the following we consider only this in-phase mode, thus 
neglecting the n dependence ^" = £,a, and = (the 
system is periodic in the x-direction with period 2Zo). We 
find that the instability appears above a flux threshold 
foc given by. 



fo > foc 



"0-1 



(9) 



where /o > 0. The stability domain in the plan (/o,<5o) 
is shown in Fig. |21 Typically the instability is related to 
a large diffusion Db on terrace b {ao 3> 1) accompanied 
by a small stiffness 7b of step Sb {So > 0). 

Although the full expression of the dispersion relation 
is cumbersome, near the instability threshold we can in- 
troduce a small parameter e measuring the distance to 
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the threshold. This parameter arises naturally when con- 
sidering the long wavelength limit, in which q —> eq. In 
this limit the relevant scaled parameters are chosen to 
be: 



/o = e/, So = eS, = I + e^a , 



(10) 



together with a rescaling of the stiffness Ta — * e^Ta. This 
scaling will lead to a consistent weak nonlinear expan- 
sion, as it will be shown below. To lowest order the two 
branches of cr = a{q, 4>) are. 



^" = -^/«(/a + 12<5)g2 - eV 



(11) 
(12) 



The growth rate ct^ is maximum for the wavenumber 
Qm = [fce{fce + 12<5)]^/^/4\/6. Even for vanishing S (no 
difference between step stiffness) an instability is still 
present driven by diffusion anisotropy and growth effect. 

Taking for example typical values for Si(OOl) at two 
different temperatures T = 1000 K and T = 800 K (with 
lo = 2 10^8 m and F = 0.1 ML s'^), /o = 0.3, and 4.3, 
ao = 45, and 119, and Sq = 0.9, we obtain that the typi- 
cal wavelength of the meanders is of the order of 250 nm 
and 60 nm, respectively. These sizes are in the range of 
the transverse modulations (ripples) of the step bunches 
observed in the experiments |20l l2lL l22l l23l | . 

We study at present the nonlinear evolution of the me- 
andering instability, in the limit of weak amplitudes and 
long wavelengths. An inspection of Eqs. (|1HI - (|12(I for 
the damping and growth rates, suggests that we should 
consider different time scales and amplitudes for the sta- 
ble and unstable branches. We introduce the functions 
s{y, t) and u{y, t) corresponding to the amplitudes of the 
stable and the unstable branches respectively. These am- 
plitudes are related to the step shape by, 



(13) 



where Mq is the matrix formed with the eigenvectors as- 
sociated to the linear dispersion relation and depends on 
the physical parameters (/o,ao:'5o)- In order to obtain 
the relevant nonlinear dynamics we use a standard multi- 
scale method. Adatoms concentrations and step shapes 
are expanded in powers of e. The amplitudes of this 
expansion depend on the slowly varying space ey and 
two time variables eH and e^t. In particular the stable 
and unstable shape functions are given by s = s{e^t, ey) 
and u = u{e^t, ey). Solving diffusion equations (|1I2|I and 
boundary conditions I|3I4|) up to order e^, and inserting 
the results into the step velocity equations H6I7|) we find 
the equation for the unstable mode: 



/ 



:^ifa + l26)u + d'yU+^{d^, 
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(14) 



where we renamed the slow variables e^t t and ey — > y. 
After rescaling we can write Eq. (|14|l in the form dtu — 
—dy\u -h dyU + {dyu)"^ /2]. This is a conserved version 
of the Kuramoto-Sivashinsky equation (CKS) which also 
describes the weak nonlinear regime of the step bunching 
instability 28] . Numerical simulations (performed with 
an extra term dyU) have revealed a non interrupted coars- 
ening dynamics, the characteristic size of coalesced step 
bunches increases as i^/^ [l^. In our context the dyU 
term, introducing a dispersive drift, is absent, but this 
should not change the i^/^ scaling. Our direct simula- 
tions confirmed this behavior and also demonstrate a lin- 
ear time growth of the characteristic meander amplitude 
{v?)^/'^ ~ t (spatial average is denoted (• • • )). A typical 
spatiotemporal evolution from a random initial condi- 
tion is shown in Fig. O In the context of Bales-Zangwill 
meandering instability, it was shown that the dynamics 
of steps is fully nonlinear, excluding CKS equation, al- 
though it would be compatible with the basic symmetries 
of the system . 

Simple similarity and matching arguments lead to a 
complete picture for the long time behavior of (|14|l . It 
is worth noting that the CKS equation admits an exact 
particular solution in the form of a stationary parabola 
u{y,t) = —y^/2. We also note that, for rapidly decreas- 
ing or bounded functions, the moment of order one of 
u{y, t) is conserved while the second order moment sat- 
isfies 



d_l 

dt 2 



u'dy^ / [(dyu)' - id'yu)']dy , (15) 



showing that the amplitude of u tends to increase in re- 
gions where the gradient term dyU inside the integral 
dominates the curvature term dyU. This suggests that 
the dynamics of large amplitude smooth regions of u{y, t) 
is almost independent of the four derivative (stabilizing) 
term in H14|l . Trying a similarity solution 



u = tXyA^), Y^y/t'^ 



(16) 



of dfU 



idyu)y2] 



we find immediately the 
exponents b = 1/2 and a = 1, which agree with our 
numerical results. The fourth order derivative term be- 
haves as dyU ~ dyf/t, which is consistent with the above 
assumption that it should be negligible in the consid- 
ered regime. Moreover, the similarity equation for ip: 
ip — (F/2)9y(y9 + dyi'P + {dy^)"^ /2\ = has a solution in 
the form of a bounded parabola: 



w(y,t) = -z/V2, \y\<yo{t) 



(17) 



and zero elsewhere. The parabola edge yo(t) is deter- 
mined by the condition 



d 1 

4, 



rvoit) 



yHy^2 y^dy, (18) 
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FIG. 3: Spacetime plot of u{y, t) with nondimensional y 
and t, given by the numerical solution of the CKS equation. 
The coarse-graining of structures leads to a superposition of 
parabolas, with a size (u^)^^^ ~ t. In the long time state all 
the parabolas tend to have unity curvature at their maximum, 
and width increasing as \/t. 



a consequence of ((T^ . which gives yo{t) = (16t/3)^/^. 
The general, asymptotic solution of H14|l can be thought 
as a superposition of parabolas satisfying H17|l (see 
Fig. |2J|. The joining region between the parabolas 
possesses a high curvature and can be described by 
the reduced inner equation dyu{y) + (l/2)(9yu(y))^ — 
2k^ = const., whose solution is of the form u{y) — 
2 log[cosh fc(j/ — 2/o(i))]. Matching with the outer solu- 
tion 1)17(1 we find that k = yo{t)/2 ^ \ft. Therefore, the 
curvature of the joining line increases as k ~ t. 

In this Letter, we have shown that the effect of step 
stiffness difference and diffusion anisotropy induces a me- 
andering instability during surface growth. We have first 
presented a linear stability analysis of a train of steps us- 
ing a simplified two-dimensional extension of the model 
studied in Ref. We have shown that the inter- 

play between diffusion anisotropy and the step stiffness 
effect under growth conditions leads to a finite wave- 
length instability which is maximum for the in-phase 
mode. Our results are complemented by a weak non- 
linear analysis of the step dynamics which reveals that 
the amplitude of the meanders is governed by the con- 
served Kuramoto-Sivashinsky equation (CKS) which dis- 
plays non-interrupted coarsening. We believe that the 
morphology observed in experiments of molecular beam 
epitaxy on Si(OOl) slightly disoriented towards the [110] 
direction, reported in Refs. [i^, 0, |U 113, can be ex- 
plained by the nonlinear evolution of the step bunching 
and step meandering instabilities arising simultaneously. 
Indeed, we will present elsewhere an investigation of the 
two-dimensional dynamics originated by the nonlinear 
coupling between these kinetic effects, and we will dis- 



cuss their role in the formation of the ripples shown for 
instance in Fig. 1 of Ref. . 
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